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S eminary 


Experiments have been carried out on the mixing of streams with and 
without swirl and with and without density differences In a coannular geo- 
metry. The flow process was visualized through the use of a phosphorescent 
tracer, ?-3 hiacetyl . Three variations of the visual Izat Ion techniques were 
used; direct excitation, collisional excitation and quenching of the tracer 
Each techrique Is complementary to the others and yields Information about 
turbulent transport, molecularly mixed and unmixed regions in the flow The 
axial variation of the flow for the case with swirl in the inner annulus is 
investigated In some detail, .‘rom the flow visualization pictures, it Is 
tentatively concluded that the Initial Helmholtz instability of the tangen- 
tial shear layer gives rise to small scale peripheral structures, then it 
Is possible that the Inner swirling stream drives an instability which turns 
these small scale peripheral structures into ones of larger scale, thereby 
rearranging the initial vorticity distribution of the flow field, The pic- 
tures also Indicate the occurrence of interchange of fluid between the 
annuli Relatively few data points were obtained for the co-flowing streams 
with density difference (nitrogen/hel ium or helium/nitrogen). During the 
course of the experimental work, an aspirating probe for concentration 
measurements was designed and built. From the flow visual ization pictures 
it is tentatively concluded that the Interaction of density field with the 
resulting radial pressure gradient may have influenced the evolution of the 


flow structures. 
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I. Introduction 

The mixing of two fluid streams, whether It be In confined or unconfined 
geometry, is of technical Interest. A unique technique of flow visualization 
had been developed at the M.I.T. Gas Turbine Laboratory by S.C. Bates *0977), 
which makes It possible to examine the molecularly mixed layer between two 
gas streams. The technique also allows one to distil ;<h the rlxed from the 
unmixed regions, S.C. Bates used the technique to study the nixing c' shear 

layer between two parallel streams. 

2 

W.K. Cheng (1978) adopted the technique to study effects of swirl 
on the mixing between two streams in a co-annular geometry. I. he course 
of his Investigation, he pointed out the existence of non-axf syrmetric 
structures In an unstable swirling flow. An effort is made in the present 
Investigation to visualize these structures In more detail for the case 
where only the Inner stream of the co-flowing streams has swirl. To achieve 

this, the optical system was modi 'led to allow better viewing, from the axial 

2 

direction, of the peripheral structures discovered by W.K. Cheng . A tunnel 
spacer was also added to enable the study of the flow field over a larger 
axial extent. An effort is also made to examine the effect of differences in 
density In such a flow system. During the course of the work, an aspirating 
probe was designed and developed for measurement of the relative concentra- 
tion of a gas of different molecular weight In a background gas, for example 
of He In Nj, this being the gas combination employed to produce the density 
difference between the annuli. The development of this technique Is In an 
early stage. 
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1 1 . Fl ow Visualization Technique 

The flow visualization technique used In the present experiment was 
developed by S.C. Bates' (1977) at the M.l.T. Gas Turbine Laboratory. 

Through the use of different modes of visualization 'the direct excitation, 
the col’ Jon, i excitation and the quenching modes), the technique can map 
ojr the ■ jlarly mixed versus the turbulently mixed regions. For details 
of the teJnlf, le, reference may be made to ref. (1). A brief description 
of the technique will be repeated here. 

The fluorescent tracer used to visualize the flow Is 2-3 biacetyl. It 
is a yellowish liquid at room temperature with a fairly high vapor pressure 

o 

(40 mm). 2-3 biacetyl absorbs light from 3500 to 4600 A, but fluoresces 
in the blue and phosphoresces in the green. The latter was used for the 
flow visualization here because it is 60 times brighter than the former. 

The quantum yield of the phosphorescence is about 15%. 

In the direct excitation mode, one of the gas streams is seeded uni- 
formly with biacetyl. The flow region of interest is illuminated by a thin 
sheet of light and the illuminated plane is photographed. Thus the image 
is a dissection of the flow at the illuminated plane. Note that unlike the 
sbadoqraph and Schlerein methods, the image is not formed as a result of 
the integrated effect on the light along the line of sight, hence its 
adaptability for visualizing three dimensional structures. Under this 
mode of visualization, a picture of the mixing region depicts the spread- 
ing of the seeded stream into the unseeded one by both turbulent transport 
and molecular diffusion. 
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In the colllslonal excitation mode, one stream Is seeded with blacetyl 
and the other with benzene. A thin sheet of UV light Is passed through the 
flow region of Interest, The UV light has no effect on blacetyl but opti- 
cally pumps the benzene. The excited benzene molecules In turn cause the 
blacetyl to phosphoresce through molecular collisions by the transfer of 
optically Induced molecular energy. Thus the recorded picture 1$ a map of 
the molecularly mixed region In the flow. 

In contrast, the quenching mode permits visualization of the molecu- 
larly unmixed region. The optical arrangement for this mode Is the same 
as for the direct excitation. One of the streams Is seeded with blacetyl 
while the other Is seeded with oxygen. In the molecularly mixed region, 
the oxygen molecules quench the phosphorescence; thus, the recorded image 
shows the molecularly unmixed region of the stream containing blacetyl. 

III. Experimental Apparatus 

A. The tunnel and the gas flow system . 

Figure 1 Is a sketch of the tunnel and the gas flow system used in the 
present experiment. A detailed description of the tunnel Is given in refer- 
ences (!) and (2). A brief description of the recent modifications is given 
here for continuity. There are two Independent manifolds Ml and M2. One 
is fed by eight cylinders of helium, and the other by twelve cylinders of 
nitrogen. These gases enter pipes A and B via the sonic orifices. These 
sonic orifices regulate the flow rates and they are so adjusted that the 
axial velocity of either stream Is about 7 ms’* at the test section (thus 
the Reynolds number based on the annulus height is of the order of lO 1 *)^. 
Ball valves VI through V4 allow the facility to be fed by gas of one type 
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only, or by gases of different types, in each pipe leading to the annulus 
tunnel, thus enabling the study of density effects. 

Gases entering the tunnel are evenly distributed over the cross-section 
by c flaking orifices on plate C, then pass through screen D into a settling 
chamber. Contraction is provided to thin out the boundary layer before 
they enter the test section. The flows (swirling if turning vanes E and F 
are mounted) emerge to form a shear layer at the end of the splitting 
cylinder G. Since oxygen quenches the phosphorescence of the experiment, 
the whole apparatus must be capable of pumping down to about 100 microns of 
mercury. It is then filled with working gas before the start of the experi- 
ment. 

The seed flow supplies from gas cylinders containing the organic liquids 
(biacetyl and benzene) diverts into the two branches. The inner annulus flow 
(A) is seeded with 2-3 bi acetyl and the outer annulus flow (B) is not seeded 
for the direct excitation experiment, seeded with air for the quenching 
experiment, and with benzene for collisional excitation experiment. The 
cylinders are heated to deliver 2 to 3 mm of partial pressure of the seed 
vapour at the test section. The long pipes A and B, with a physical length 
to diameter ratio of about 120 ensure that the organic chemicals are uni- 
formly mixed before entering the tunnel. 

Because the tracer gas 2-3 bi acetyl used in the experiment has a very 
unpleasant odor, the flow has to be discharged via a pneumatically controlled 
plexiglass door (also serves to keep air out) into a huge plastic baloon. 

A blower is then used to blow the tracer contaminated gas through an elec- 
tric kiln (fig. 2) which has been heated to a temperature of about 2000°F. 
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so that the organic chemicals may be broken down before discharging into 
the building exhaust system. Five Instrument ports are provided at the 
tunnel test section to allow the taking of hot-wire data and aspirating 
probe data. They are two Inches apart with the first one 3 inches down- 
stream from the splitting cylinder. 

A test section spacer Is provided so that the structure of the shear 
layer may be examined over a larger axial distance. 

B. Input Output Optics 

Details of the input and output systems are in references 1 and 2. 

A 200 Joule capacitive flashlamp pulser drives an EG and G FX-38-C-3 
Xenon flash tube to provide the Input sheet of light. The light Is colli- 
mated by a slit and two cylindrical quartz relay lenses Into a thin sheet 
before entering the tunnel through a UV grade quartz window. A Pomfret 
Research Optics 2537 A Interference filter Is used to obtain the UV light 
for the collislonal experiment. The transmission efficiency of the filter 
Is about 20% and the bandwidth is about 280°A. Figure 3 shows the optical 
arrangement. 

The dim light emitted in the visualization experiment, after being 
collected by a Nikon f = 55mm/ 1:1. 2 lens, passes through a green filter and 
falls on the Input area of an ITT F 4747 mlcroc.iannel image intensifier. 
Newly purchased, the gain of this device was 2x10"*. Over the years, with 
continuous use, its gain has deteriorated. The resulting ouput Image of 
the flow region is recorded through contact photography with Polaroid 
type 57 film {ASA 3000). Pictures for direct excitation and quenching 
experiments are taken with exposure times of 750 us at f/2.0, while those 
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for colllslonal mode are taken with exposure time of at least 1.5 ms at f/1.2. 

C) Ouratlon of the experiment . 

The tunnel works In a blow down mode to conserve gas. The flow is 
initiated by a pneumatically actuated valve from a reservoir of gas cylinders. 
It lasted for about 2.5 sec. From the time history of the pressure tracer 
(Fig. 4) at the mixing pipes and tunnel, the flow is steady for at least 1.5 
sec. 


D) H ot wire measu re ments . 

A triple hot-wire probe (DISA type 55-P91), operating In constant tem- 
perature mode. Is used for measuring the three components of flow velocity at 
the test-section. Time resolved data, obtained by radially traversing the 
probe, or with the probe at a stationary point, are acquired over an axial 
distance of 18" (with the use of test section spacer). The duration of data 
acquisition is about 0.7 sec., during which the flow is steady. The data 
acquisition system is shown in fig. 5. The signal from the hot-wire is 
digitized using an MP 6912 analog-to-digital converter (with a capability 
of sampling at 100 KHz). A user built logic circuit then transfers the 
digital signals to the core memory of a DEC MINC 11/2 microcomputer, via a 
parallel input and output board on the microcomputer Q-bus. About 16,000 
points of data are collected in each run. Thus a sampling frequency of 7 
KHz is used for each hot-wire. After each run, the data are transferred 
from the core memory to a secondary storage (floppy diskettes) for per- 
manent storage. 
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Previous documentation (Including the computer program) of the use of 
this hot-wire probe was lost, therefore they had to be redeveloped for the 
present experiment. Some details for the use of this hot-wire probe are 
presented In Appendix A. 

E. Aspirating probe measurement . 

An aspirating probe has been designed to measure the relative concen- 
trations of the constituents of binary mixture of gas. Details of the de- 
sign of this probe and its operation are in Appendix B. The probe is esti- 
mated, using axisynmetric analysis, to have a spatial resolution of 0.65 inn 
and a frequency response of 25 KHz . When tested in a shock tube, the output 
signal had a rise time of 40 us, as shown in the Appendix. The probe was 
also tested in a stream of nitrogen of varying velocity for its response 
to changes in external velocity. This test indicated no changes in the out- 
put readings of the anemometer, verifying its design concept. 

Some results are reported from this probe, but it must be considered in 
an early stage of development. 

IV. Experimental data. 

A, Introduction 

In this section some experimental data are presented for the case 
where the inner annular stream has swirl, while the outer has none. The 
gas used for streams of equal density is nitrogen, while those for streams 
of different densities are nitrogen and helium. For each run, a flow 
visualization picture and a set of hot-wire probe traverse data/ or an as- 
pirating probe traverse data can be obtained. Though efforts are made to 
ensure reproducibility in conditions for runs with the same nominal values. 
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variations in them may occur. Relatively few data points were obtained 
(24 pictures and 4 sets of raw data from the aspirating probes) for the 
co-flowing streams with different densities, all with the inner one swirling. 
In contrast to reference 2, flow visualization pictures of axial planes are 
taken with the camera mounted in such a way that it looks into the flow 
axially (see fig. 3b for the optical arrngement). It was hoped that this 
would enable better viewing of the peripheral structures of the flow-field. 

A t jnnel spacer was added to the system to enable the viewing of the flow 
field over an axial distance of 14 in. Also it must be mentioned that, in 
contrast to the flow visualization data presented in references 1 and 2, 
only the direct excitation and quenching mode pictures are presented in 
details. The collisional excitation pictures obtained were of poor quality. 
At this stage, it is tentatively concluded that this is largely due to the 
deterioration of the gain of the ITT F4747 microchannel image intensifier 
through its usage over the years. 

B. Hot-wire data. 

Hot-wire traverse measurements were made at 10 axial stations. The 
stations are two inches apart (except 5 and 6, which are 1" apart), with the 
first one three inches downstream of the splitting cylinder. Fig. 6 shows 
the mean velocity components (radial, tangential and axial) plotted as func- 
tions of radius at axial stations 3", 5", 7", 9", 11" downstream of the 
splitting cylinder. These mean velocity components were obtained from 
280 point averages (a simple moving average) over 5120 data points. There- 
fore these represent a time domain average of 4 ms and space domain average 
of 5 mm. It should be noted that the accuracy of the velocity components as 
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measured by the triple wire probe is subject to the limitations of the method 
used, as described in Appendix A. Furthermore, the triple wire probe is cali- 
brated for speed and direction under static conditions but used for measure- 
ment of velocity components in highly turbulent flow. 

No attempts are made to retrieve the fluctuating velocity components 
from the hot-wire data. We present some of the linearized signals from the 
triple wire probe in Fig. 7, for stations 1 to 10. These signals contain 
information aboj„ the flow structures. All the signals shown in Fig. 7 rep- 
resent a fuV' ' . pie-wire probe traverse outwards, starting at 1 cm 

from the centre body. At station 1, the signal has higher fluctuations in 
the inner annulus than in the outer annulus, where the signal is compara- 
tively quiet. Note the discontinuity in the signal {in the neighbourhood 
of mid-annulus), at stations 1 and 2. This discontinuity indicates the 
position of the shear layer. The large fluctuations in the inner annulus 
seem to be spreading outwards progressively. At stations 9 and 10, the hot- 
wire signals seem to indicate that the flow is turbulent throughout the 
annulus. Careful perusal of the hot-wire data indicates that there is a 
definite sawtcoth envelope superimposed on the fluctuating signal; this is 
especially well-defined at station 5. This probably indicates the presence 
of some recirculating eddies. 

C. Data from aspirating probe. 

Aspirating probe traverse measurements were made at axial stations 
1 and 4 for two flow configurations, (i) where the inner annular flow was 
nitrogen and the outer annular flow was helium, and (ii) where the flow 
configuration in (i) was reversed. The plots in Ficures 8 and 9 show 
the mean density profile for the four traverse measurements. The output 
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shown is the anemometer readings in volt, rather than in concentration. 

This is because the hot-wire element within the probe was broken while 
making the measurements; consequently, we are unable to check the calibra- 
tion. 


D) Flow visualization pictures . 

Flow visualization pictures are presented for two cases; 

(i) the co-flowing streams of equal density with the inner one set to 
swirl at a flow angle of 40 to 50°, (ii) the co-flowing streams of unequal 
densities with the inner one set to swirl at same flow angle as in case 
(i). For case (ii), nitrogen is used in one stream while helium is used in 
the other. Some flow visualization data were obtained for the cases with 
the inner stream set to swirl, while the outer stream was set to swirl in 
the opposite direction and with the swirl in both streams, but they are not 
presented here. Pictures of the flow are taken on Polaroid type 57 film. 

The circular image is limited in area by the 18 mm diameter image intensifier 
output face (but they can be blown up). All the flow visualization pictures 
presented here are dissections of the flow field in a radial plane parallel 
to the axial direction or in an axial plane perpendicular to the axial direc- 
tion. The center of the picture corresponds approximately to the middle of 
the annulus. For the pictures in a radial plane, the flow is from the right 
to the left. 

The degree to which information can be retrieved depends very much on 
the quality of the flow visualization pictures. Reference may be made to (1) 
for discussion of this topic. The gain of the image-intensif ier is spatially 
ncn-uniform, the image being considerably dimmer near the edges, This factor 
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should be taken into consideration while attempting to interpret the data. 

The collisional excitation pictures obtained are of poor quality; this is 
tentatively attributed to deterioration In the gain of the image intensi- 
fier unit. 

Case ( 1 ) : 

In this case, both the Inner and outer annular streams were of nitrogen. 
The inner stream was set to swirl at an angle of 45 ° to 50°. The axial 
velocity was about 7 to 8 ms"^ . Flow visualization pictures were taken 
over an axial extent of 14". No extensive or intensive efforts are made 
here to Interpret each and everyone of the pictures obtained here, but fea- 
tures of the flow field that can be adduced frcm the pictures will be indi- 
cated. 

The pictures numbering from 1 to 6 in fig. 10 are the direct excitation 
pictures of the flow field in the radial pla'ie parallel to the axial direc- 
tion. The number in inches by the side is the distance in inches the center 
of the picture is from the splitting cyli ne'er. The letters a, b and c refer 
to flow visualization pictures taken at the same location with the flash 
tube fired at different times after the flow is steady for the blow-down 
period. This procedure is followed with the tacit assumptions that (i) the 
flow Is repeatable from run to run with the experimental conditions kept 
constant (evidence may be provided by ? ome of the flow visualization pictures), 
(ii) there may exist some structural eddies which are convected with the 
flow (note that the flow visualization picture represents a three inch cir- 
cular flow field). 
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Pictures (Cl a and #1 b have similar overall global features. The undu- 
lating Interface may Indicate the beginning of Helmholtz Instability. How- 
ever, pictures #2a to 02c, 03a to 03c, show some striking features. Picture 
02a seems to Indicate motion of fluid in the inner stream outwards (perhaps 
due to the centrifugal action of swirl, in accordance with Rayleigh's criter- 
ion) ^hlle pictures 02b, 02c, 03a and 03c, indicate the entrainment of lumps 
(very loose technical word) of fluid from the outer stream into the inner 
stream. What causes such transport of fluid elements and what role does 
such mechanism play in the mixing of the two streams? This could be due to: 
(i) the results of Hemholtz instabil ity ,or (ii) the results of the swirling 
flow which rearranges the vorticity present in the inner swirling flow into 
large scale motion, as seen in these pictures. Such features are not pre- 
sent in the flow visual ization pictures for linear shear layers in references 
(1) and (2). Therefore it is te ntativ ely concluded that the resulting in- 
stability (as given by Rayleigh criterion) may have resulted in the rear- 
rangement of the vorticity originally in the inner swirling stream into such 
eddy motion. Another picture worth noting is number 5a. This picture indi- 
cates that some of the fluid elements have been pulled away entirely from 
the inner stream. Picture 05b indicates the fluid from the inner stream 
may have permeated the outer stream. Quenching pictures (fig. 11), indi- 
cate that at such a distance from the splitting cylinder, molecular mixing 
is rather complete among fluid elements from both streams. Pictures 01 to 
05 in fig, 12 are the collisional excitation pictures. They show such poor 
definition that no attempt will be made to discuss them. 
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Pictures 1 to 15 In fig. 13 (direct excitation mode), show the dissec- 
tion of the flow field In the axial plane at different distances from the 
splitting cylinder. Picture #1 shows a rather quiescent Interface between 
the streams, while pictures 2 and 3 Indicate some non-axl symmetric flow 
structures (some kind of small eddies) have begun to form. There could be 
a large number of such azimuthal structures around the peripheral. It Is 
possible that the formation of such a large number of azimuthal modal struc- 
tures may precede the development of large azimuthal modal structures, such 
as are evident In pictures 6, 7 and 10. Such large azimuthal modal structures 
are still present In pictures 11a, 12, 13, and 14b. It Is of interest to 
note that some fluid elements from the outer stream have penetrated into 
the inner one near the centre body, as is evident in pictures 10 and 11 
This could indicate an exchange of fluid elements between the two streams 
has occurred. Pictures 12, 13 and 14 are indicative of the fluid motion 
in pictures 2 and 3 of fig. 11. Pictures 14 and 15 indicate that fluid 
elements from the inner stream have permeated the outer stream and some 
fluid elements have actually parted from the inner stream. Pictures #1 
to #10 in fig. 14 are the corresponding quenching pictures. Picture #10 
of fig 14 and the collisional picture #1 of fig. 15 indicate that molecu- 
lar mixing of eddies may have been quite thorough at 9-1/2" from splitting 
cylinder. From this set of pictures, we may tentatively conclude that the 
swirl in the Inner stream (giving rise to instability, in accordance with 
Rayleigh's criterion) may have been responsible for the large scale motions 
In the neighborhood of the interface between the two streams, and in par- 
ticular for the non-axisyrrmetric structures. The small scale peripheral 



15 


structures formed close to the splitting cylinder may be the result of 
Helmholtz Instability of the tangential shear layer caused by swirl. The 
centrifugal action of the Inner stream could then drive an Instability, 
giving rise to larger peripheral structures. 

Case (11) - 

Some preliminary flow visualization pictures have been obtained for 
the co-flowing streams of different densities (one helium and the other 
nitrogen) with the inner one set to swirl at 45° to 50°. The helium stream 
is set to have an axial velocity of about 8 to 9 ms"^ at the test section, 
while the nitrogen stream has an axial velocity of about 6 ms - ^ at the test 
section. For this case, visualization pictures were taken only over an 
axial extent of 7". These are shown In figures 16 to 23. The pictures do 
indicate some qualitative differences in the flow features for the two 
configurations: (a) nitrogen in the inner stream and helium in the outer 

stream, (b) helium In the inner stream and nitrogen in the outer stream. 

For the co-flowing streams of different densities with the inner one 
set to swirl, the following may play a part in the evolution of the flow 
field: (i) the swirl which sets up a pseudogravity field, (ii) the gravity 
effect set up by the unequal density interface (the strength of such effect, 
though small, varies circumferentially), and (iii) the difference in the fluid 
properties of the two streams (e.g., viscosity) which may influence the micro- 
structure of the flow field. The influence of (i) and ( i 1 ) may reinforce one 
another (e.g., at the lower half of the tunnel test section) or counteract 
against one another (e.g., in the upper half of the tunnel test section with 
nitrogen in the outer stream and helium in the inner stream.) 
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In the following, we will make an attempt to Interpret the data for the 
two flow configurations, bearing In mind that there may be too few data 
points to come to any conclusive evidence of any particular flow features 
unique to such a flow system. 

(a) with nitrogen In the Inner stream and helium In the outer stream: 

We have noticed that picture 1 of fig. 10 exhibits a continuous undu- 
lating feature (possibly Indicating the Inception of Helmholtz Instability 
as mentioned), while such Is not evident in picture 1 of fig. 16, except for 
the protruding fluid element into the outer stream. 

It is possible that the density difference at the interface under the 
influence of a radial pressure gradient In the inner stream, could have 
accounted for the observed difference In flow features between the two. 

Local stability analysis of the shear layer with a radial pressure gradient 
in the inner stream shows that such a flow configuration would give rise to 
an enhanced Rayleigh-Taylor instability; this is indeed evident in pictures 
2 and 3. In picture 3, the nitrogen has spread into the outer annulus while 
helium from the outer annulus has penetrated far into the Inner annulus. 

Fig. 17 shows the corresponding quenching pictures. Picture 3 shows an 
almost unmixed strand of fluid element lying across the annulus. Such a 
feature is not evident in all other quenching pictures obtained so far. 

What is seen ma^ be a result of fluid motion driven by a density gradient 
in a pseudo-gravity field (i.e., the presence of a pressure gradient). The 
vortlcity equation for an inviscld fluid with variable density is 

ar + V • Vfi = (£1 . V) V + ~ (Vp x 7p) 


( 1 ) 
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where .1 is the vorticity, V the velocity, p the density and p the pressure. 
This equation says that In the presence of pressure gradient (or acceleration 
field), vorticity Is rr*»*ted in a fluid of variable density In motion. For 
such flows, vorticity, or rotation, Is created so as to rotate the fluid 
until surfaces of constant density are normal to the pressure gradient (or 
acceleration vector field). Thus the resulting fluid motion might have 
rotated the fluid element from the Inner annulus Into the position shown In 
picture 3. It must be pointed out that the above explanation is arrived 
at through speculation on the mechanisms which influence evolution of fluid 
flow in streams of equal and different density and by process of elimination. 

The pictures In figs. 18 and 19 are the direct excitation and quench- 
ing pictures taken for flow field in the axial planes. Pictures 1 and 2 
in fig. 18 shows ft rmatlon of small eddy-like azimuthal structures, while 
picture 3 (though overexposed) indicates the formation of a larger azimuthal 
structure. As for the quenching pictures, picture 3 shows that a consider- 
able fraction of the fluid in the inner stream is molecularly mixed out 
(contrast with picture 4 of fig. 14 for homogeneous flow). 

(b) with helium in the inner stream and nitrogen in the outer stream: 

Picture 1 of fig. 20 indicates that the interface is not as quiescent 
as that in picture 1 of figs. 10 and 16. Further downstream, swirling 
instability (in the sense of Rayleigh's criterion) develops, thus enabling 
the spreading of the fluid from the inner stream to the outer stream as is 
evident In pictures 2 and 3 of fig. 20. The pictures in fig. 21 are the 
corresponding quenching ones. Notice the strand of almost unmixed fluid 
element from the inner stream lying across the annulus (though not as distinct 
as that in picture 3 of fig. 17). As in the case of configuration in ii, this 
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Is likely to be due to the fluid motion driven by the density gradient In 
a combined gravity-acceleration field. ) 

The pictures In figure 22 and 23 are the direct excitation and quench- 
ing pictures of the flow field In the axial plane. They again Indicate the 
development and presence of non-axl symmetric peripheral structures. 
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V. Some Concluding Remarks 

Attempts were made to Investigate the flow field of a co-flowing annu- 
lar flow system over an axial distance of 14" through flow visualization 
pictures and hot-wire data. The Inner annular stream had swirl, while the 
outer one did not. The effect of density difference on such a flow system 
Is touched upon. An aspirating probe has been designed and built for concen- 
tration measurement, but no extensive use has been made of it for measurements 
yet; as only four traverse measurements of the concentration in binary gas 
flow system (helium/nitrogen) have been made. Interpreting the resulting data 
involves some extrapolation and speculation. It is tentatively concluded 
that: 

(i) the initial Helmholtz instability of the shear layer gives rise to 
small scale peripheral structures; it is possible then that the inner swirl- 
ing stream drives an instability (in the sense of Rayleigh's criterion), 
which turns these small scale peripheral structures into ones ot larger 
scale, thus rearranging the initial vorticity distribution of the flow, 

(ii) the hot-wire data indicates that the spread ot turbulence into 
the outer annular stream is progressive rather than abrupt, 

(iii) some of the hot-wire data have well defined envelopes, possibly 
indicating the presence of recirculating eddies. 

It is also tentatively concluded that the density difference has an 
effect on the overall structure in such a flow system. Basing on the rela- 
tively few data points we have obtained so far, it is tentatively concluded 
that for the co-flowing steams of equal and unequal density in the present 
experimental set up, any qualitative difference in the flow features 
observed further downstream of the splitting cylinders is more likely caused 
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by the generation of vortlcity due to the density gradient In the presence 
of a pseudo-gravity field (i.e., the radial pressure gradient). 

VI. Suggestions for future work . 

In the present set up for the optical system, one can only get a 2 to 
3 inch circular view of the flow field. Based on such data, one then attempts 
to make a guess of the several flow structures. The optical arrangement 
should be modified to permit viewing of the larger flow field, especially 
a larger peripheral extent from the axial direction. This will aid in 
clarifying the flow structures. 

More data have to be obtained for the co-flowing streams with unequal 
densities so that the effects of density difference can be elucidated. 

The aspirating probe can be used in conjunction with the triple wire probe 
to measure the velocity components. 

There are many other parameters, such as those mentioned In reference 
2; for instance, the Reynold number, which might affect the flow structures. 
They have yet to be Investigated. 

It is noted that the annular flow system Is of simple geometry and thus 
it makes the spectral method (i.e., pseudospectral , Tau or Galerkln Approxi- 
mations) attractive for numerically simulating the flow structures in such 
a system. Such method for studying the ' asic fluid phenomena in such a 
flow system may be considered in conjunction with future experimental work. 

Lastly, future efforts should also be directed in an attempt to quantify 
the experimental data. 
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Appendix A 

The Triple Hot-wire Probe Operation 


(a) Introduction 

The use of the DISA type 55-P91 triple hot-wire probe goes along the 
same line as described In Appendix A of reference (2) with minor difference 
in the method of retrieving velocity components. However, the computer 
program developed In reference (2) for fitting calibration data to models 
representing the operation of the triple probe and those for retrieving 
the three components of velocity from the hot-wire readings cannot be 
traced, consequently new computer programs for the purposes have to be re- 
developed for the present experiments. The probe is operated in constant 
temperature mode using the three identical anemometers described in A2 of 
Appendix A in reference (2). In using this probe, it is assumed that the 
dependence of speed and geometry of the triple wire probe can be neglected; 
otherwise its use would involve a thorough calibration of the probe at ell 
possible orientations and within speed range of interest at each orientation. 
This assumption also permits the user to do the speed calibration and 
directional calibration separately. 

( b ) C al ibratio n 

The calibration (both speed and directional) of the triple wire probe 
was done in the low turbulence wind tunnel of the M.I.T. Aero & Astro 
Department. The M1NC 11/2 microcomputer was used for the purpose. The 
readings from the triple hot-wire and that of pressure transducer (dynamic 
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head of wind tunnel section) were digitized using the Datel ST-LSI-2 
analog-to-dlgltal converter. At the end of calibration, they were trans- 
ferred to floppy diskettes for permanent storage. 

(1) Speed Calibration 

Each wire in the probe is calibrated against the speed in the wind 

tunnel by orientating the probe such that the output reading from the 

( 2 ) 

anemometer Is the maximum . This ensures that the wind direction is per- 
pendicular to the wire. The wire is calibrated in air but the result is re- 

(21 

scaled for use in nitrogen ’ using King's Law: 

Nu = A + B r£ (A.l ) 

or 

hr = A + B (^) n ( A - 2 ) 

where Nu = Nusselt No 

Re - Reynold No 

E = Output voltage from anemometer 

K = thermal conductivity of gas 

p = density of gas 

p = viscosity of gas 

The constant A and B are found from the calibration data using linear re- 
gression analysis. 

Each wire was calibrated for speeds ranging from 0 to 15 ms"^. 


at the mean temperature 
of the wire and the 
ambient 
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(11) Directional Calibration 

Let be the effective wind speed as seen by wire i (1=1 to 3) 
and u = (uj, u 2 , u 3 ), where are the velocity components In a suitable 
laboratory co-ordinates system, then the effective velocities are re- 
lated to u through 


where 

l 


UJ 

Ht 2 

l 2 2 

Z 3 2 


U3 [«] + [ cl 


and Q 


ui 2 " 
Uj 2 

u 3 2 


UiU 2 


U 2 U 3 ! 
UjUl 


The and C are the calibration constants given by 


A lt 

Am 

Aj 3 

Ai 4 

Ai 5 

Aj 6 ' 


‘ c, 

A2 1 

A2 2 

A 2 3 

A2 

A2 5 

A26 

and C = 

C 2 

As 1 

A3 2 

A 3 3 

A34 

A 35 

A 3 6 


■ c 3 . 


(A. 3) 


Thus altogether there are twenty-one calibration constants. They 
are found by fitting the model to the m calibration data points using 
multiple linear regression analysis, which yields a matrix equation. They 
can be inverted to yield A and C, i.e., 

f\t % 



24 




mV] 

'Xi 


"1 


[ * ] T [ 6 t ] 


(A. 4) 


where 



[ <*, ] ■ 

[ A 1i A i 2 A i 

, \ A i s 

\ 

c i3 


l*i 1 * 

[(y)i . . . 

( Z 1 2 )j ’ * 

■ <v>»] 

and 

‘ (U; 2 ), 
• 

« 

<U 2 2 )x (u 3 Z )l 

(UlU 2 )! 

(u 2 u 3 )i 

(UlUi)* 

[ x j - 

<- ! >0 

4 

(u 2 2 )j ( U 3 2 ) j 

(UiU 2 )j 

(u 2 u 3 )j 

( u 3 U ! ) j 



(u, 2 ) (u 3 2 ) 

' 2 'm v 3 ; m 

(u lU2 ) m 

< U!U =>m 

(u, Ul ) m 

where i runs from 1 to 3 

i, L 

and j denotes the j tn calibration points 




The directional calibration is done by orientating the probe for a few 
sets of speeds. A typical example would involve rolling the triple-wire 
probe through 240° for 3 sets of yaw angles and 2 to 3 sets of speeds (in 
the range of interest), thus m can run from 126 to 225. 
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(111) Retrieving the velocity components: 

The equations (A. 3) relating the effective velocities as seen 
by the triple wire probe to the velocity components In the laboratory co- 
ordinates form a set of non-linear algebraic equations. Given the three 
values of effective wind speeds (as determined from the triple hot-wire 
pro..' readings), we have to solve the set of non-linear algebraic equations 
for the three-components of velocity in laboratory co-ordinates system. 

It is the author's experience that sometimes as many as eight sets of roots 
may be obtainable. Thus, the flow direction has to be partially known for 
retrieving the velocity components from the knowledge of effective wind 
speeds. Here, we use the Levenberg/Marquardt iteration method (modification 
of classical Newton-Raphson iteration) to look for the three components of 
velocity from eqn. (A. 3) (in contrast to the first-order gradient search 
method in reference 2). 

The method goes as follows: 


(k) 

We replace ' estimate ot u by 


• v (k) * ^ k) 

( k) 

where Au is the solution of the system of linear equations 


(A. 5) 


3 


l 

j=1 


(Y 


(k) T 


ij 


3 

t=l 


(k) 

ti 


J (k) l 
°tj } 


Au 


(k) 



(A. 6) 
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where k is the iteration index, I is the identity matrix and J is the 

n, * 

Jacobian. 


- [?] 


1 


*\t 


(k) 

>Jw k) 



U 






.00 


The Jacobian j}|0 cai either be evaluated numerically or analytically, 
is a positive parame<er which is being calculated so that 



( “ k+1) ] 2 < 


5 t l( u< k >)] 

i=l 1 * 


(iv) Some remarks, n the use of triple hot-wire proble: 

t 

Some results of tJTe directional calibration of the probe are shown in 
figures (A.l) to A.4)£ As pointed out in Appendix A of reference (2), the 
model equation (A.3l Jjoes give negative values of and has been forced 

4 

to zero in the figures. Also there are regions where the calibration 

curves are co-1 inear, which may cause an uncertainty in the determination 

/ 2) 

of the velocity components. 

Solutions u to equation (A. 3) are non-unique. It seems that to use 

r \j 

this triple probe for measuring the fluctuating components of velocity, one 

may have to resort to a thorough calibration of the probe. The probe also 

suffers from drift in cold resistance and thus it has to be calibrated 

(2) 

every time before use. ' a time consuming process. 

It must also be mentioned that in using the Levenberg/Marquardt itera- 
tion scheme, the method may converge to a stationary point which is not a 
solution. 
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Appendix B 
Aspirating Probe 

(i) Introduction 

An aspirating probe was designed and built with the aim of measuring 
gas concentrations in co-flowing streams of differing density. In what 
follows, a brief description of the operating principle and Its design is 
given. 

(ii) Operating Principle 

Basically, the probe (Fig. B.l) consists of a hot-wire sensing element 
(tungsten is used here) placed in a cylindrical duct of constant area fol- 
lowed by a throat. The throat exits into a chamber being pumped to vacuum, 
so that with the probe placed In a stream, the flow Is choked at all times. 
The hot-wire sensing element is connected by leads to an anemometer { DISA 
55001) operating in a constant temperature mode. The mass flow rate is 
given by 


(pU> 5 ■ 



(rui 

ir Pt 

2 2(y-1) 

<9§T> 

- 

" JV 


where p is the density 
U Is the velocity 

P is the pressure 

T is the temperature 

Y is the specific heats ratio 


n 4 , i _ universal gas constant _ Ro ^ 

R Is the gas constant ( S pecu la r weight > 

and subscripts t and s refer to stagnation values and sonic values, 
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Thus for a given gas mixture, the maximum flow rate per unit area is depen- 

p 

dent upon t^r only, ^ the stagnation conditions are held constants, 
then for a given gas, the Mach number M is fixed for a fixed geometry, 

(1 .e. , fixed area ratio, 

Using the basic gas dynamic relations, it can be deduced that p, p, T and 
U are constants, so that the hot-wire element responds only to gas concentra 
tlon and stagnation conditions. For a binary mixture of gases, at the same 
T^, the Pj, can change due to the viscosity of gases (especially at solid sur- 
faces, or mixing). But the fractional change in can at most be 0(M J ). 

Thus for low speed flow, such a probe may only respond to gas concentration. 

In the next section, we will indicate how one can estimate the response 
time of the probe. 

(Hi) Estimation of Response Time 

In this analysis, it will be assumed that 

(a) the lateral velocity v is much smaller than the longitudinal velocity 
U H.e., V < < „ 

(b) the flow is incompressible up to the probe intake 

(c) the flow is inviscld 

The method of singularities (sinks, sources and ring vortices) allows 
the flow past this probe to be calculated. Instead of using distribution 
of ring vorticies (a more elaborate and complicated way involving the mani- 
pulation of various types of integral), wc use a much simpler approach here. 
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The probe may be approximated by a half-Rankine axisymmetrlc body. The 
sampling action of the probe may be simulated by a sink olaced at the tip 
of the probe (See fig. B,2). If the sampling rate of the probe is much 
smaller than U», CU 2 , then the Ranklne axlsymmetric body will not be altered 
significantly. With these In mind, and with reference to flg.B.2, the stream 
function for the flow Is given by 

V X {x ' a) 


X = ir z U» + lL-g 




- u« a 2 






(0.2a) 


where d^ Is the diameter of the capture streamtube, and a = ^ D™, with D* 
as the diameter of the unperturbed Rankine-type axlsymmetric body. In dimen- 
sionless form (based on 0» and U,,}, the streamlines are given by 


, <*- v 
!. r 2 _ + J_ 

2 16^ + (x-1) 7 ' 8 




+ r 


7 1 


= constant (B.3) 


while the stagnation points are given by (r,x) which satisfies 


1 + 


(* - ?) 


V 


16 [r 2 + (x- p 2 ]h 8 ^ x2 + r2 ^ 


= 0 , 


(B.4) 


and 


16 


V 


l 1 mx - 1 y] H 8 0(2 + r2)>! 


• 0 


(B.fi) 


simultaneously. 
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The above axlsymmetrlc analysis would give a flow field as shown In 
fig. B.3, For the dimension of the probe shown In fig. B.l, the lateral 
spatial resolution of the sampling volume would be about 0.7 mm. The lon- 
gitudinal spatial resolution of the sampling volume Is determined by the 
frequency response of the probe (e.g., for a flow speed of 7 ms"^ with a 
probe response of 25 KHz, the resolution would be about 0.30 rrm). Fig. 

B.4 shows the variation of time t taken by the fluid element in the sam- 
pling volume to reach the probe with the streamlines within the capture 
stream tube. Of course, it would take the fluid element on the stagnation 
streamlines an Infinite amount of time to reach the probe. Thus, an esti- 
mate of the response time t would be given by 



f* 


For the given built probe shown in fig. B.l and Fig. B.5, the response 
time was estimated to be about 40 ps. The given probe was tested in the 
M.I.T. Gas Turbine Laboratory shock tube, as shown in Fig. B.6. The rise 
time was measured (using Nicolet Explorer Digital Oscilloscope Model 206) 
to be about 40 ps (Fig. B.7a and B.7b). 

It is noted that the sensing element is placed at a position in the 


duct where the flow (laminar) is expected to be fully-developed. 
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fig. 2; Electric kiln- use# to break down tie 
organic chemicals before dumping into 
the laboratory exhaust system. 








Upper beam: at mixing pipe of 

inner annulus 

Scale: Horizontal = 0.5 sec/cm 
Vertical = 0.2v/cm 

CO psig/cm) 

Lower beam: at test section 

Scale: Horizontal - 0.5 sec/cm 
Vertical = 0.05v/cm 

(0.5 psid/v) 


Fig. 4a: Flow-time history (with no tunnel spacer added on} 



Upper beam: at mixing pipe of 

outer annulus 

Scale: Horizontal = 0.5 sec/cm 
Vertical = 0.05v/cm 

(2.5psig/cm) 

Lower beam: at test section 

Scale: Horizontal =0.5 sec/cm 
Vertical = 0.05v/cm 

(0.05psid/cm) 


Fig. 4b: Flow-time history with tunnel spacer added on. 
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I Analog signal from hot-wire/ 
t aspirating probe and position indicator 



DEC MI INC 11/2 microcomputer 


Fig. 5 Data acquisition system 
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Fig. 6a Mean velocity components at station 1 
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Fig. 6b Mean velocity components at station 2 
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Fig. 6c Mean velocity components at station 3 
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Fig. 6d Mean velocity components at station 4 
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Fig. 6e Mean velocity components at station 5 


otLoaiTY in li/t 

O. 1. 


««m vctositT 



%• *••<•** 

jtAoiua in sn 



44 


Fig. 7: Linearised signal from the triple-wire probe at 

station 1 through station 10 
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Fig. 8a 
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Fig. 8b 
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Fig. 9a 
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The pictures are blown-up view of the originals using Polaroid 
type #52 film. 


Notation: 

The letter following the picture number (e.g., #la) refers 
to flow visualisation picture taken at the same location 
with the flash-tube fired at a different time after the 
flow is steady for the blown-down period. 

The number in seconds is the time at which the flash-tube 
s fired. 


The number in inches is the distance the center of the 
picture is from the splitting cylinder. 



53 


Fig. 10: Direct excitation pictures (as viewed from side window) 

for co-flowing streams of equal density with the inner 
stream set to swi rl . 

Direction of flow is from right to left. 


Data Format 


Outer annulus 



Approximate initial position of shear 
(also of density discontinuity if present) 
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Fig. 13: Oirect excitation pictures of flow field in axial plane 

(as viewed from an axial direction) for co-flowing 
streams of equal density with the inner stream set to 
swirl . 


Data format 
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suction rate * 1.59 x 10~ m 3 /s 
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free stream velocity * 7 ms" 



Fig. 6.3: Flow around an ax i symmetric Rankine-type body with a 

weak sink placed at its tip. 









Fig. B. 7a : Output of Digital Qscillo? 

response time of probe “ t 
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